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By introducing a predictive mechanism with small-world connections, we propose a new motion 
protocol for self-driven flocks. The small-world connections are implemented by randomly adding 
long-range interactions from the leader to a few distant agents, namely pseudo-leaders. The leader 
can directly affect the pseudo-leaders, thereby influencing all the other agents through them effi- 
ciently. Moreover, these pseudo-leaders are able to predict the leader's motion several steps ahead 
and use this information in decision making towards coherent flocking with more stable formation. 
It is shown that drastic improvement can be achieved in terms of both the consensus performance 
and the communication cost. From the industrial engineering point of view, the current protocol 
allows for a significant improvement in the cohesion and rigidity of the formation at a fairly low 
cost of adding a few long-range links embedded with predictive capabilities. Significantly, this work 
uncovers an important feature of flocks that predictive capability and long-range links can compen- 
sate for the insufficiency of each other. These conclusions are valid for both the attractive/repulsive 
swarm model and the Vicsek model. 

PACS numbers: 05.65.+b, 89.75.-k, 89.20.Kk 



I. INTRODUCTION 

Over the last decade, physicists have been looking for 
common, possibly universal, features of the collective 
behaviors of animals, bacteria, cells, molecular motors, 
as well as driven granular objects. The collective mo- 
tion of a group of autonomous agents (or particles) is 
currently a subject of intensive research that has po- 
tential applications in biology, physics and engineering. 
One of the most remarkable characteristics of systems, 
such as flocks of birds, schools of fish, and swarms of 
locusts, is the emergence of collective states in which 
the agents move in the same direction, i.e., an ordered 
state [H, H, i, i, i, S 0). Moreover, this ordered state 
seeking problem for flocks/swarms/schools can be fur- 
ther generalized to consensus [8j, rendezvous, synchrony, 
cooperation and so on. From the application aspect, this 
kind of distributed collective dynamic systems has direct 
implications on sensor network data fusion, load balanc- 
ing, unmanned air vehicles (UAVs), attitude alignment 
of satellite clusters, congestion control of communication 
networks, multi-agent formation control and global coor- 
dination for emergency d, Qj], [HI, G3 ■ 

The interaction pattern of the natural biological 
flocks/swarms are neither entirely regular nor entirely 
random. An individual of a flock usually knows its neigh- 
bors, but its circle of acquaintances may not be confined 
to those who live right next door. In 1998, in order to 
describe the transition from a regular lattice to a random 
graph, Watts and Strogatz (WS) introduced the concept 
of the small-world network [131 ] by rewiring one end of 
a few connections to new nodes chosen at random from 
the whole network. With these few shortcuts, the aver- 
age distance is decreased significantly without crucially 
changing the clustering property. The work on the WS 



small- world network has started an avalanche of research 
on complex networks, especially, the synchronizability 
of networks can be greatly enhanced by introducing a 
few long-range connections [3, EH [IB]- Thus, for bet- 
ter synchronization in a flock of neighboring-connected 
agents with a leader, it is advantageous to build a small- 
world-type network structure by randomly adding long- 
range connections from the leader to a few distant agents 
(namely pseudo-leaders), so that the leader can affect 
them, thereby influencing all the other agents through 
them, via fast communication and rapid control com- 
mands. 

Although a lot of relevant works were focused on net- 
work structures, recently, more and more researchers are 
interested in finding the rules of the inter-connections 
present in abundant bio-groups. Extraction of these rules 
can help interpret why the bio-groups can demonstrate so 
many good characteristics such as synchronization, sta- 
bilisation, cohesion, etc. A fairly basic but popular flock- 
ing strategy can be traced back to the Reynolds Model 
171 ] , in which three elementary flocking protocols are pre- 
scribed, (i) separation: steer to avoid crowding and col- 
lision; (ii) alignment: steer towards the average heading; 
(iii) cohesion: steer to move towards the average posi- 
tion. These rules have been proven to be effective and 
thus become the basic rules for the design of bio-group 
dynamic models. In 2003, Gazi and Passino [l8| proposed 
an effective A/R (attractive/repulsive) swarm model in 
which the motion of each individual (autonomous agent 
or biological creature) is determined by two factors: (i) 
attraction to the other individuals on long distances; (ii) 
repulsion from the other individuals on short distances. 

The Gazi and Passino A/R model [Tq |. embedded with 
a similar mechanism of the inter-molecular force, is de- 
rived from the biological flocks/swarms behaviors. Thus 
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far, the general understanding is that the swarming be- 
haviors result from an interplay between a long-range 
attraction and a short-range repulsion between the in- 
dividuals In [ijj], Breder suggested a simple 
model composed of a constant attraction term and a re- 
pulsion term which is inversely proportional to the square 
of the distance between two members, whereas in [20j 
Warburton and Lazarus studied the effects on cohesion 
of a family of attraction/repulsion functions. Moreover, 
in physics communiuty, a large volume of literature on 
systems with interactive particles have also adopted func- 
tions of attractive and repulsive forces to investigate the 
dynamics of the system [H 11 H HI SE HE [p] . For 
instance, in [2lL [22l | , systems of particles interacting in a 
lattice are considered with attraction between particles 
located at different sites and repulsion between particles 
occupying the same site. It is discussed in [23j ] that, the 
structure of a nonuniform Lcnnard-Jones (LJ) liquid near 
a hard wall is approximated by that of a reference fluid 
with repulsive intermolecular forces in a sclf-consistently 
determined external mean field incorporating the effects 
of attractive forces. 

With the A/R model [18j| . Gazi and Passino proved 
that the individuals will form a bounded cohesive swarm 
in a finite time. One year later, by adding another fac- 
tor, i.e., attraction to the more favorable regions (or re- 
pulsion from the unfavorable regions), they generalized 
their former model into a social foraging swarm model 
[28| . Under some suitable circumstance, agents in this 
modified model are apt to move to the more favorable 
regions. The A/R model has been adopted by physicists 
and biologists to model self-driven particles and biolog- 
ical flocks d [H HE SH HI Hi In 2004, Moreau |j 
presented a linearized model of flocks, and proved that 
the flock is uniformly and globally cohesive to a bounded 
circle if and only if there exists an agent (the "leader" ) 
connecting to all other agents, directly or indirectly, over 
an arbitrary time interval. Other than these kinds of 
A/R models [H, HE [34]], another popular kind of model 
is that without leaders (say homogeneous), where a very 
representative one is the Vicsek Model Q. In each step, 
every agent updates its velocity according to the average 
direction of its neighbors. With the decrease of external 
noise or the increase of the particle density, the collective 
behavior of the flock undergoes a phase transition from 
the randomly distributed phase to the coherently moving 
phase. In 2003, Jadbabaie et al. provided the conver- 
gence condition of the noise- free Vicsek model, i.e., the 
individuals are linked in some intervals 
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Although most of the previous works on flock dynam- 
ics yield many advantages such as synchronization, sta- 
bilization, cohesion, and quick consensus, agents within 
the networks only know the information that is currently 
available to them. In this paper, we highlight another ap- 
pealing phenomenon, i.e. the universal existence of pre- 
dictive mechanism in various biological aggregated sys- 
tems. A general physical picture behind this is illustrated 
in Fig.[T]and interpreted as follows: in widely-spread nat- 
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FIG. 1: (Color online) Predictive vision in natural bio-groups. 
Each individual makes its motion decision based on not only 
the current status of the leader and its neighbors but also 
their future dynamics. 



ural bio-groups composed of animals, bacteria, cells, etc., 
the decision on the next-step behavior of each individ- 
ual is not solely based on the currently available state 
information (including position, velocity, etc.) of other 
(neighboring) agents inside the group but also on the 
predictions of future states. More precisely, taking a few 
past states of its leader and neighbors into account, an 
individual can estimate the corresponding future states 
several steps ahead and then make a decision. 

Some experimental evidences have already been re- 
ported in the literature. In 1959, Woods implemented 
some experiments on bee swarms and found a certain 
predictive mechanism of electronic signals inside this bio- 
group [lU. Also for bee swarms, in 2002, Montague et 
al. discovered that there exist some predictive protocols 
in the foraging process in uncertain environments [37| . 
Apart from the investigation of the predictive mecha- 
nisms of locust swarming and foraging, more scholars fo- 
cused on the predictive function of the optical and acous- 
tical a ppa ratuses of the individuals inside bio-groups 
[38l[39l.T4Cl| . especially cortexes and retinae. For instance, 
based on intensive experiments on the bio-eyesight sys- 
tems, they found that when an individual observer pre- 
pared to follow a displacement of the stimulus with the 
eyes, visual form adaptation was transferred from current 
fixation to the future gaze position. These investigations 
strongly support our conjecture of the existence of some 
predictive mechanisms inside abundant bio-groups. 

Bearing in mind the plentiful examples of predic- 
tive protocols inside natural bio-groups, we incorporated 
some predictive functions into a few long-range links, and 
found that it is possible to significantly enhance the flock- 
ing performances at a fairly low cost of the additional pre- 
dictive energy. More interestingly, proper prediction ca- 
pability can help reduce the minimal number of the long- 
range links between the leader and the pseudo-leaders, 
thus effectively decrease the communication cost. 

On the other hand, from the industrial application 
point of view, the phenomena and strategy reported in 
this paper may be applicable in some relevant prevail- 
ing engineering areas like autonomous robot formations, 
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sensor networks and UAVs 0, E3, E3l- Typically, due 
to the limitation of the communication energy, only a 
few agents have the capability to communicate with the 
leader. The incorporation of a predictive mechanism 
into these pseudo-leaders can greatly improve the flock- 
ing performances. 

The rest of this paper is organized as follows. In Sec- 
tion HH the small-world connection model with embed- 
ded predictive mechanism is presented. Then, in Sec- 
tion IIII1 its important role in improving flocking syn- 
chronization performances is extracted and analyzed by 
numerical simulations on the A/R model [HI]. After- 
wards, in Section IIV1 the generality of the virtues en- 
dowed by suchpredictive mechanism is validated on the 
Vicsek model [3| as well. Finally, conclusions are drawn 
in Section [V] 

II. MODEL 
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FIG. 2: Small-world predictive mechanism of flocks. Here, the 
leader (L), pseudo-leaders (P), and followers (F) are denoted 
by star, triangles and circles, respectively, and the neighboring 
area of each individual is a circle with radius r centering at 
itself. The leader's trajectory is given in advance, which will 
not influenced by the others. Each pseudo-leader's dynamics 
is always influenced by both the leader's future position H p 
steps ahead and its neighbors' positions, while each follower's 
dynamics is solely affected by its neighbors' positions. 



It is well-known that a ring-shaped network structure 
is not a good one for efficient mutual communication and 
global control within a flock of agents, while the so-called 
small-world networking structure performs much better. 
By adding a few long-range connections, the average path 
length of the ring-shaped network will be abruptly de- 
creased. This small- world effect is very desirable for 
fast communication and information transmission, effi- 
cient synchronization, and effective global control over 
the entire network [151 ] . Thus, in a flock of neighboring- 
connected agents with a leader, for communication and 
control purposes, it is advantageous to build a small- 
world-type network by randomly adding long-range con- 
nections from the leader to a few distant agents (namely 
pseudo-leaders) . As shown in Fig. [2 the leader can affect 
pseudo-leaders, thereby influencing all the other agents 
through them. To be clear, we call the non-special agents 
followers. Thus in our model, there are three different 
kinds of agents: leader (L), pseudo- leaders (P), and fol- 
lowers (F). 

In this model, the flock is assumed to move in an Tri- 
dimensional space and the standard A/R function [§.[l8l. 

m 

G{d pL ) = -d pL (a - b ■ exp(-||dp L ||!/c)) (1) 

is used as long-range interaction from the leader (L) to 
each pseudo- leader (P), where a, b, c are three free pa- 
rameters, d P L is the m-dimensional vector pointing from 
the predicted location of leader L to the current location 

of a pseudo-leader p, and ||<ipz,||2 = J dpjjdpL denotes 

the Euclidean distance between them. The force G(d P L) 
is an m-dimensional vector whose direction is from the 
pseudo-leader to the leader. For simplicity, in our model, 
the motion of leader is given in advance, and will not be 
affected by any other agents. We assume every pseudo- 
leader has the same prediction horizon H p , that is to say, 
a pseudo-leader will predict the leader's location H p steps 
ahead. 




FIG. 3: (Color online) Information communication process 
inside flocks with a predictive mechanism. The arrows repre- 
sent the passing of position information. If the pseudo-leader 
i predicts the dynamics of the leader H p steps ahead, then 
the follower j, within topological distance equaling H p from 
i, could be affected by the current location of the leader. Here, 
an agent j has topological distance zero to itself, and 1 to all 
the agents located inside the circle with radius r centered on 
j. Two agents having topological distance 1 are called con- 
nected. An agent has topological distance Dt to j if and only 
if it is connected to at least one agent with distance Dt — 1 
to j, while is not connected with any agents having distance 
smaller than Dt — 1 to j. Note that all the agents are not nec- 
essarily embedded in a circular pathway. Instead, due to the 
circle-shaped neighborhood of each agent, the position infor- 
mation transmission pathway forms a tube tangent to these 
neighboring circles, and the layout boundary of the pathway 
could be any types of curves depending on the location of the 
agents. 

On the other hand, a weaker A/R function, repre- 
senting the short-range interaction between two arbitrary 
neighboring agents i and j is addressed as: 

g(dij) = -dij fa-&-exp(-||dy||§/c)) , (2) 

where dij is the m-dimensional vector pointing from the 
individuals j to j, and ||<iij||2 = \J dfj denotes the Eu- 
clidean distance between them. The parameters a, b and 
c are much smaller than a, 6, and c, respectively. The 
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direction of vector g(dij) is from i to j. Denote by r the 
radius of neighboring area (see Fig. \2§ ■ The neighboring 
A/R links could connect any two agents (F-F, P-P and L- 
F) within the Euclidean distance r except the L-P inter- 
action described in Eq. (fTJ). Note that the leader can in- 
fluence other agents, but will not be influenced. In order 
to decrease the prediction cost, no predictive mechanism 
is incorporated into the neighboring A/R links. Bearing 
in mind the physical meaning of A/R function the 
positions of a pseudo-leader z p and a follower Zi (both z p 
and Zi are m-dimensional vectors) are determined by 

z p (t) = G(d pL (t + H p )) + Yl 9(d pj (t)), (3) 

lonq link to the leader j^^-dpj (t)_ r 

neighboring links 

and 

m= E 9(<iijm (4) 

j,dij(t)<r 

\ ■* 

neighboring links 

respectively, where t denotes the current time, and 
d P L(t + H p ) represents the m-dimensional vector pointing 
from the leader's position H p steps ahead to the current 
position of a pseudo-leader. Although the rest of this pa- 
per concentrates on the motions in 2-dimensional space, 
the present model can be directly applied in any finite 
dimensional space. In this way, unlike the routine flock- 
ing strategies [H, [H, [H, |4l[ , a small- world interaction 
pattern is established with embedded predictive mecha- 
nism, which has the capability of predicting the future 
behavior (position, velocity, etc.) of the leader several 
steps ahead. Note that the structure of this interaction 
network will change in time since the location of each 
agent is varying. 

The information communication process is illustrated 
in Fig. [3J The farthest agent i% directly communicat- 
ing with agent i is among the ones at the rim of the 
circle with radius r centered on agent i. Analogously, 
the farthest agent ii directly influenced by i\ is also lo- 
cated at the rim of the circle centered on agent i±, and 
so forth. Finally, the influence of agent i reaches agent 
j in H p steps. When agent j receives the information 
from agent i at time t, it is in fact a delayed informa- 
tion of agent i at time step t — H p . However, if agent i 
acts as a pseudo-leader who can accurately predict the 
behavior of the leader H p steps ahead, then, at time step 
t, agent j's motion is affected by the exact current loca- 
tion of the leader zi,(t). In this way, although agent j 
may not have direct connection with the leader, it could 
know some information of the leader's current dynam- 
ics by agent i's delayed information. Therefore, agent 
j can adhere to the leader more tightly, the flock's for- 
mation is more likely to be stable, and the coherence of 
the whole flock is thus improved effectively. Note that, 
the predictive mechanism is valid only if the leader's mo- 
tion is regular [42j. If the leader moves in some random, 



chaotic or other irregular ways, such as random walk, it 
is, in principle, impossible for the pseudo-leaders to pre- 
dict the leader's further location. Fortunately, in the real 
biological world, the flock leader always moves in some 
predictable pattern. Therefore, the other agents have 
the opportunity to display their predicting ability, such 
as that used by a chameleon to capture a fly and by a 
dog to catch a frisbee. 

III. ANALYSIS AND SIMULATIONS 

To show the advantages of the predictive mechanism, 
we compare the performances of the two cases of flocking 
with and without the predictive mechanism by simula- 
tions over an A-agent flock moving in a two-dimensional 
space as shown in Fig. [J] The parameters are set as fol- 
lows: the neighboring circle r — 0.65, the parameters of 
the A/R functions (see Eqs. fTJ) and ^j) of long-range 
links and neighboring links are set as a = 8, b = 17.6, 
c = 3.2, and a = 1, b = 2.2, c = 0.2, respectively. The 
former A/R function is much stronger in order to inten- 
sify the influence of the leader. As shown in Fig. [JJi, 
each agent starts from a position randomly selected in 
the square [0,1] x [0,1]. The leader and the pseudo- 
leaders are selected randomly among these N agents. 
The trajectory of the leader is set as xi = y/xi, and 
the velocity of the leader is vl xi (t) = 0.02, Vi, (t) = 

y/Q.02(t + 1) — V0.02i. In our simulations, the time la- 
bel t is a discrete number with step length being equal 
to 1. 

It can be seen from Figs. HJHIJI that, with a proper 
H p , the coherence of the flock will be improved remark- 
ably. The followers will adhere to the leader much more 
tightly (see Fig.[4j:), and the flock formation will be more 
stable. More precisely, for the flocks with the predictive 
mechanism, the position error index 

1 N 

i=\.i^L 

converges to a constant after finite steps, indicating a 
stable state of the flock dynamics. Here, J p measures the 
cohesion performance of the flock, with 1 1 rfi_c 1 1 2 denoting 
the Euclidean distance between agent i (F or P) and the 
leader. Meanwhile, as to the flock without this mecha- 
nism, as shown in Fig. [4|d, J p will keep increasing along 
with the elapse of time, making the flocking unstable. 
However, abusing the foresight, namely over-prediction 
(see Fig. 2J1), is also undesirable. That is because the 
pseudo-leaders are attracted/repelled by the leader posi- 
tion too many steps ahead and will probably escape the 
flock with a fairly high speed, and then lose influences on 
the followers. In this way, the flocking will be damaged 
after finite steps. 

The circular formation of the followers in Figs. Hfc-Hli 
is because of the particular form of the A/R function. To 
be clear, we give a more detailed explanation as follows: 
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FIG. 4: (Color online) (a) Starting position of the 50-agent 
flock consisting of 1 leader, 2 pseudo-leaders and 47 followers. 
The blue line marked by square points denotes the trajectory 
of the leader, (b) Flock position after 65 steps without predic- 
tive mechanism (H p = 0), (c) with proper predictive mech- 
anism (Hp = 20), and (d) with over prediction (H p = 70). 
Here, x\ and X2 denote the two-dimensional position coordi- 
nates. 
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FIG. 5: (Color online) The roles of the pseudo-leaders' num- 
ber N p i (figures (a) and (c)) and prediction horizon H p (fig- 
ures (b) and (d)) on a flock of 50 agents. The leader and 
the pseudo-leaders are selected randomly among these agents. 
Each point is an average over 1000 independent runs. The pa- 
rameters of the A/R functions ((TJ) and © are a — 8, b — 17.6, 
c = 0.4, and 5=1,6= 2.2, c = 0.2, respectively. The radius 
of the influence circle is r = 0.65. Each agent starts from a 
position randomly selected in the square [0, 1] x [0, 1]. With- 
out loss of generality, the trajectory of the leader is set along 
the curve defined by X2 = y/xi, and the velocity of the leader 



Actually, the biological flock dynamics is fairly complex. 
For example, fish uses sidetrack to sense the current vari- 
ances, in this way, the fish schools are formed. On the 
other hand, in order to save flying energy for the rear 
individuals, wild geese flocks always form a 'A'-like for- 
mation [43(1 . In this case, the forming process of such flock 
is determined by aerodynamics. Our proposed model is 
based on the idea that an individual inside a bio-group 
can predict the trajectory of its leader (s), and this kind 
of intelligence can help the individual make more efficient 
decision to improve the flocking performance. This pa- 
per, however, does not aim at reproducing the detailed 
movement formations of any particular bio-groups. 

In order to extract the role of H p and the number of 
pseudo-leaders denoted by N p i, we display their influ- 
ences on the position error index J p and velocity error 
index J v in Fig. 5, where 
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Here, J v measures the formation performance of the 
flock, where vt and vt denote the velocity vectors of the 
leader and the ith agent (F or P). If J v — » 0, the relative 
velocity of each pair of agents approaches zero, thus the 
flock formation is fixed. In Fig. [5^,, we fix H p and dis- 
play the curves of J v with increasing N ph while Fig.[5jD, 
on the contrary, reports the curves of J v with increas- 
ing Hp and fixed N p i . It can be seen from Fig. [5]a, that 
the curves fall sharply at the beginning and then more 
slowly until reaching a minimum, afterwards rising slowly 
as the increasing of N p i. It implies that adding just very 
few pseudo-leaders (e.g. long-range links) to the leader, 
which transforms the flock topology from a strongly lo- 
calized network into a small-world one, will improve the 
flocking performance greatly. However, when the number 
of pseudo-leaders reaches an optimum N* t corresponding 
to the minimal J* , the flock formation performance will 
start to worsen and these extra pseudo-leaders become 
redundant. On the other hand, increasing H p can help 
reduce J v in two ways: (i) it depresses J v with the same 
N p i; (ii) it reduces the optimal value of A^*; correspond- 
ing to the minimal J* . Compared with Fig. [5^i, the J v 
curves in Fig. [5]b fall more slowly at the beginning until 
reaching the minimum, afterwards J v increases all along 
and never reaching a stable platform. It implies that 
the flock formation performance can be remarkably im- 
proved with proper predictive capability, however, too 
much vision into the future, namely over-prediction, will 
even worsen the formation of flocking. On the other 
hand, more pseudo-leaders, i.e., larger N p i (as long as 
Npi < N*{), can also help yield more cohesive flocking 
with better formation. 

An interesting phenomenon can be observed from 
Fig. [5ji that increasing the number of pseudo-leaders be- 
yond some characteristic number makes the flocking per- 
formance worse, which seems counter-intuitive. This can 
be explained as follows. First, as shown in Eqs. ^ and 
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FIG. 6: (Color online) Illustration of Fig. [5ji. Follower C 
heads in the red dashed arrow direction with the guide of 
pseudo-leader A who moves nearly in the same direction as 
the leader, making the movement direction nearly synchro- 
nized. However, the addition of pseudo-leader B will make 
the direction of C deviate from the red dashed arrow to the 
green arrow, and thus worsen the velocity synchronization 
performance J v . Therefore, more pseudo-leaders are not nec- 
essarily beneficial to synchronizing the velocities of the flock. 



insight into the future and moderate number of pseudo- 
leaders are preferred. 

A plausible physical rule behind the observed phenom- 
ena shown in Fig. [5] is that, in order to achieve a fixed 
flocking performance, greater predictive capability and 
more pseudo-leaders can compensate the insufficiency of 
each other. Given a fixed-size flock, separately for forma- 
tion performance J v and cohesive performance J p , there 
exists an optimized combination of H p and N p i . The con- 
clusions achieved in this paper are not sensitive to the 
trajectory of the leader. In order to validate the gener- 
ality of the conclusions, we have used another two differ- 
ent trajectories, i.e., parabolic and sinusoidal curves. The 
simulation results are shown in Appendix B, which sug- 
gest that our main conclusion, i.e., predictive capability 
and long-range links can compensate for the insufficiency 
of each other, also holds in those two cases. 



(TjJ , the velocities of both the pseudo-leaders and the fol- 
lowers are solely determined by the relative positions. 
Indeed, as shown in Fig. [6l if there is only one pseudo- 
leader A, then the follower C will head in the red dashed 
arrow towards A, and all the three individuals are nearly 
moving in the same direction, which is desirable. How- 
ever, when another pseudo-leader B is added into this 
flock, the follower C will move towards a certain position 
between A and B (see the green arrow) . Apparently, the 
velocity synchronization performance is worsened com- 
pared with the former case. Therefore, the observation 
in Fig. [5^ is reasonable. To give a more vivid and detailed 
explanation, we give further discussions in Appendix A 
for interested readers. 

Next, we investigate the effects of H p and N p i for an- 
other important index J p . It can be seen from Fig. [5fc 
that the curves fall sharply at the beginning and then 
asymptotically approach a stable value. The main dif- 
ference between Fig. [Sh. and Fig. [5J; is that the latter is 
monotonous and has no minimum, in other words, the 
increase of N p i always improves J p . A special case ex- 
plaining this phenomenon is that if all the followers serve 
as the pseudo-leaders, then they will be very cohesive to 
the leader. However, when Npi exceeds a certain value 
N p i , J p will increase so slowly that almost no substantial 
improvement can be achieved. Moreover, increasing H p 
can help reduce J p , and the improvement shrinks along 
with increasing H p , indicating a saturation effect that no 
remarkable improvement can be achieved by pushing H p 
to a very high value (which may cost too much). Com- 
pared with Fig. [St, the curves of Fig. decrease more 
slowly first and then reach the lowest values at a fairly 
large H p , afterwards rise quickly. Thus, over-prediction is 
not preferred. Analogous to Fig.[S]b, more pseudo- leaders 
will help improve the cohesive flocking performance, how- 
ever, this improvement also displays a saturation effect 
with increasing N p i. Actually, when N p i exceeds a cer- 
tain value, J p increases very slowly that almost no benefit 
can be gained by further increasing N p i . In brief, suitable 
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(a) (b) 

FIG. 7: (Color online) With-noise case vs. noise-free case for 
a flock with N = 50, N pt = 9, and r/ = 0.02. The other 
parameters and the initial conditions are the same as those 
in Fig. [5] Each point is an average over 1000 independent 
runs. The moderate external noise does not change the global 
behavior of the flock. 

From the practical point of view, external perturba- 
tions (noise) and internal modeling mismatch are always 
present in any realistic systems, which inevitably induces 
some prediction error for the pseudo-leaders. To ex- 
amine the influence of such prediction errors, we now 
introduce the perturbation into Eqs. (TjQ) and @ by 
adding an external two-dimensional white noise term 
£ € [— \rj, ^rj\ x [—\rj, ^rj] to the vector d p L, i.e., d p L — 

d P L and G(d pL ) = -d pL (a - b ■ exp (— 1| d pi II 1 / c ) ) ■ 
The leader's position is no longer perfectly known to the 
pseudo-leaders. From Fig. one can observe that mod- 
erate external noises do not change the global behavior 
of the flock. The tendency of the curves J p and J v are 
almost the same as the noise-free case. Furthermore, to 
understand the capacity and robustness of our purposed 
predictive mechanism more deeply, we have also investi- 
gated the influences of stronger prediction errors £ with 
other kinds of leader trajectories including parabolic and 
sinusoidal curves in Appendix C. We found that the 
tolerance range of prediction error is large enough. In 
this way, the generality of the conclusions on the role of 
the predictive mechanism is thus further verified. 
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IV. PREDICTIVE MECHANISMS IN THE 
VICSEK MODEL 

The role of predictive mechanisms highlighted in Sec- 
tion [HT] is not merely confined to A/R flocks but quite 
general. To verify this, we now incorporate this predic- 
tive mechanism into another flocking model, i.e., the Vic- 
sek model [3|, and compare the synchronization perfor- 
mance of the predictive small-world Vicsek model with 
the one of the standard Vicsek model. 




(a) 




W /N 
pi 



FIG. 8: (Color online) (a) Snapshot of the predictive Vicsek 
flock at the 12 th running step. The red particle denotes the 
leader; the green particles represent the pseudo-leaders and 
the blue particles denote followers. The centered black circle 
outlines the trajectory of the leader. In this case, the predic- 
tion horizon is H p = 4. (b) Velocity synchronization index J v 
as a function of the parameters H p and N p i . The parameters 
of this simulation are L = 15, 77 = 0.1, v — 0.15, N = 300, 
r = 1, and R = L/6. The simulation result is averaged over 
1000 independent runs. 

In this model, the velocities Vi of the N agents com- 
posing the group are determined simultaneously at each 
discrete-time instant, and the position of the ith agent is 
updated according to 

x t (k + 1 ) = a;,- (ft) + Vi(k), 

where Vi (ft) denotes the velocity vector of agent i at time 



ft. For each agent the velocity vector, Vi(k), is character- 
ized by a constant magnitude v and by a direction 9i(k) 
whose dynamics is given by 

e i {k + l) = {6 i {k)) r + A6i, 

where (9i (k)) r denotes the average direction of all the 
agents' velocity vectors within a circle of radius r cen- 
tered on agent i, i.e., 



(*)>r = 



arctan [(sin (0< (ft))) r / (cos {9 l (ft)))J 

if (cos(0i(k))) p > 0; 
arctan [(sin {9 l (ft))),, / (cos (6>; (fc))) r ] + tt 

otherwise, 



where (sin (9i (ft))) r and (cos (0j (ft))) r denote the aver- 
age sine and cosine values, and A9i represents a ran- 
dom noise obeying a uniform distribution in the interval 
H7/W2]. 

As shown in Fig. [8^ the particles are distributed in 
a square of dimension [0, L] x [0, L\, The trajectory of 
the leader, which is not affected by others, is a circle 
centered at (L/2,L/2) with radius R = L/6 so that the 
direction of the leader changes constantly. The small- 
world predictive connection framework shown in Fig. [5] 
is used together with the Vicsek model. Hence, the N p i 
pseudo-leaders are always influenced by the leader's ve- 
locity Hp steps ahead together with its neighbors' current 
velocities. It is shown in Fig. [5Jd that drastic improve- 
ment of the velocity synchronization performance can be 
achieved with moderate prediction horizons. Similar to 
the results of the A/R model shown in Section 11111 onc 
can also conclude that suitable insight into the future 
and moderate number of pseudo-leaders is preferable. 



V. CONCLUSION AND DISCUSSION 

Inspired by the predictive mechanisms that universally 
exist in abundant natural bio-groups, we incorporate a 
certain predictive protocol into flocks with small-world 
structure. The predictive mechanism embedded in the 
pseudo-leaders lets many neighboring-connected follow- 
ers know the current or even future dynamics of the 
leader in time, thus each individual can make a deci- 
sion based on timely information of the leader instead 
of the delayed information as in some traditional mod- 
els. In this way, the followers become more cohesive to 
the leader and the flock formation becomes more stable. 
Note that, in our model, the leader's motion governs the 
trajectory of the whole swarm. However, it is a gen- 
eral feature of real migration flocks. Actually, in [4l|, 
a changeable target known by a few leaders is used to 
guide the whole flock, in which only the target's motion 
drives swarming behavior. Analogously, in this paper, a 
single leader is used to determine the general trajectory 
of the whole swarm. Therefore, the current leader-driven 
swarm is not unrealistic. 
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Simulation results led to the following conclusions: (i) 
Increasing the number of pseudo-leaders can always im- 
prove the cohesive flocking performance. Furthermore, 
it can improve the formation flocking performance when 
the pseudo-leader number has not exceeded a threshold, 
otherwise, the performance will be degraded, (ii) There 
exists a certain value of H p that optimizes the cohesive 
and the formation flocking performances, in other words, 
moderately increasing H p will improve the flocking per- 
formance, whereas predicting too many steps ahead will 
impair the flocking, (iii) Predictive capability and long- 
range links can compensate for the insufficiency of each 
other. It is worthwhile to emphasize the observed over- 
prediction phenomenon, which is of significance in prac- 
tice. 

Furthermore, to verify the generality of these conclu- 
sions, we have also applied the predictive mechanism to 
another popular flock model, the directed graph model 
with linear dynamics [44l . l45| . The corresponding results 
also strongly suggest that predictive protocols are benefi- 
cial to flocking dynamics when taking into consideration 
both the flocking performance and the communication 
cost. More importantly, with this mechanism, only a 
very small proportion of the followers are required to act 
as the pseudo-leaders to achieve a better flocking per- 
formance, as measured by J p and J v . From the indus- 
trial application point of view, the value of this work 
is two-fold: (i) The flocking performance is significantly 
improved by injection of a suitable predictive mechanism 
into the pseudo- leaders; (ii) Moderately increasing the 
predictive capability can help remarkably decrease the 
required number of pseudo-leaders. The latter feature 
is fairly useful for networks with insufficient long-range 
communication links, which are routinely costly. 

This work provides a starting point aimed at achieving 
better flocking performance by using a predictive mech- 
anism, and we hope that it will open new avenues in this 
fascinating direction. 
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Appendix A: Supplemental illustration of Fig. [5b 

We take a 4-agent flock for instance. As shown in 
Fig. [5^,, all the individuals except the leader are pseudo- 
leaders. However, due to the repulsion between each cou- 
ple of pseudo-leaders and the powerful attraction from 




(b) 



FIG. 9: Supplemental illustration of Fig. [5^. In this scenario, 
pseudo-leader A moves nearly in the same direction of the 
leader, and degrading pseudo-leader C (see sub- figure (a)) 
into a follower (see sub-figure (b)) helps improve the velocity 
synchronization performance J v . 

the leader, the pseudo-leaders will keep certain distances 
between each other and form a section-like shape. Since 
all the pseudo-leaders always point to the leader's predic- 
tive position Hp steps ahead, the moving directions of all 
the individuals are obviously far from synchronized. By 
contrast, when one pseudo-leader C has been degraded 
into a follower as shown in Fig. [5Jd, and take into consid- 
eration that the followers are always lagging behind the 
pseudo-leaders, the direction of follower C is aligned to 
approach the heading of the leader more or less. In this 
way, the velocity synchronization performance J v is im- 
proved. Thus, too many pseudo-leaders are not preferred. 
For more general cases with prediction errors, more long- 
range connections will introduce larger deviation, which 
partially neutralize the positive effect of pseudo-leaders 
on J v , thus moderate pseudo- leaders are also desirable. 

Appendix B: Effects of the predictive mechanism 
with different leader's trajectories 

The conclusions drawn in this paper are not sensitive 
to the trajectory of the leader. In order to validate the 
generality of the conclusions on the role of predictive 
mechanisms, we have used another two different trajec- 
tories, i.e., sinusoidal (x2 = sin(2xi) + 1) and parabolic 
(x2 = x\) curves. As shown in Fig. 9 and Fig. 10, 
no matter what the leader's trajectory is, the influences 
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FIG. 10: (Color online) The roles ol the pseudo-leaders' num- 
ber N p i (figures (a) and (c)) and prediction horizon H v (fig- 
ures (b) and (d)) on a flock with a total of N = 50 agents. 
The trajectory of the leader is set along the sinusoidal curve 
defined by X2 = sin(2a;i) + 1, and the velocity of the leader is 
(*) = 0.02, vl^ (t) = sin(0.04(t + 1)) - sin(0.04t). The 
other parameters and initial conditions are the same as those 
in Fig. [5] 





(d) 



FIG. 11: (Color online) The roles of the pseudo-leaders' 
number N p i (figures (a) and (c)) and prediction horizon H p 
(figures (b) and (d)) on a flock with a total of N = 50 
agents. The trajectory of the leader is set along the parabolic 
curve defined by X2 — x\, and the velocity of the leader is 
v Lxi (t) = 0.02, v L „ 2 (t) = (0.02(£ + l)) 2 -(0.02i) 2 . The other 
parameters and initial conditions are the same as those in 

Fig.ua 



of Hp and N p i on the principal tendencies of both J p 
and J v remain the same. Those simulations suggest that 
our main conclusion, i.e., predictive capability and long- 
range links can compensate for the insufficiency of each 
other, is also valid in those two cases. 



Appendix C: Effects of prediction errors 
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FIG. 12: (Color online) The effects of prediction error £. Here, 
the trajectory of the leader is set along the curve defined by 
X2 = y/xl (sub-figures (a) and (b)), X2 = sin(2xi) + 1 (sub- 
figures (c) and (d)) and X2 = x\ (sub-figures (e) and (f)), 
respectively. The flock size is N = 50, and the number of 
pseudo-leaders is N p i = 9. The other parameters and the ini- 
tial conditions are the same as those in Fig. 5. Each point is 
averaged over 1000 independent runs. It can be found that 
moderate prediction error £ does not change the global behav- 
ior of the flock. If £, however, reaches a very large value, like 
r\ — 0.5, the benefits of prediction capability on the synchro- 
nization J v will almost vanish. Note that, the curve of J v is 
more sensitive than that of J p , making the feasible prediction 
error range of J„ smaller than the counterpart of J p . 

In realistic system models, mismatch and external 
perturbations which cause prediction errors are always 
present. For completeness, we hereby examine the ef- 
fects of the prediction errors on synchronization be- 
havior of the flock. As shown in Fig. [T^l we present 
curves of J v and J p with increasing noise magnitude r\ = 
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0, 0.02, 0.2, 0.5 for the cases of three different leader tra- 
jectories (square root, parabolic and sinusoidal curves), 
respectively. It is observed that moderate prediction er- 
ror £ (£ < 0.2, for example) can hardly change the syn- 
chronization tendency of the flock, which guarantees the 
feasibility and superiority of the current predictive mech- 
anism. On the contrary, too large £ (£ > 0.5, for exam- 
ple) will inevitably impair the advantages of this predic- 
tive mechanism. This observation is reasonable, since 



every method has its own limits and one can not expect 
a poor prediction capability to yield a superb guidance 
for the flock. Fortunately, the tolerance range of predic- 
tion error is satisfactorily large, which further verifies the 
generality of our proposed predictive mechanism. Note 
that the curve of J p is more robust than the curve of J v , 
which roots in that the relative velocities can be more 
easily deviated by prediction error than the relative po- 
sitions. 
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